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We report e"'"e~ — > bb cross section measurements by the BABAR experiment performed during an 
energy scan in the range of 10.54 to 11.20 GeV at the PEP-II e'^e~ colMder. A total relative error 

of about 5% is reached in more than three hundred center-of-mass energy steps, separated by about 
5MeV. These measurements can be used to derive precise information on the parameters of the 
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T(10860) and T(11020) resonances. In particular we show that their widths may be smaller than 
previously measured. 

PACS numbers: 13.25.Hw, 14.40.Nd 



Recent discoveries of non-baryonic charmonium states 
that do not behave as two-quark states [l[ call for a search 
for other resonances befonging to this possible new spec- 
troscopy. Given the charmonium content of these new 
states, one could infer the presence of similar resonances 
containing b quark pairs. The observed J^'-' = 1 ex- 
otic states (r(4260), y(4350), and y(4660) 0) scaled up 
by the mass difference between the J/-!/" and the T{1S) 
(AM ~ 6360 MeV/c^) would be exotic bottomonium 
states with masses above the T(4S') and below 11.2 GeV. 
Moreover, the T(10860) and the T(11020) states, which 
are candidate T{5S) and T(6S) respectively, were ob- 
served in the same region [3, [i| . 

Between March 28 and April 7, 2008 the PEP-II e+e" 
collider Q delivered colliding beams at a center-of-mass 
energy {^/s) in the range of 10.54 to 11.20 GeV. First, 
an energy scan over the whole range in 5 MeV steps, 
collecting approximately 25pb~^ per step for a total of 
about 3.3 fb^^, was performed. It was then followed by 
a 600pb~^ scan in the range of -y/s=10.96 to 11.10 GeV, 
in 8 steps with non-regular energy spacing, performed 
in order to investigate the ^(§S) region. This data set 
outclasses the previous scans [3, |j] by a factor > 30 in the 
luminosity and ~ 4 in the size of the energy steps. Across 
the scan, the energy of the positron beam was kept fixed 
at 3.12 GeV, while the electron beam energy was varied 
accordingly, to set the required ^/s. This produced a 
variation of the boost of the center-of-mass frame during 
the scan. 

In this Letter we present, for each step in ^/s, the mea- 
surement of Rb{s) = crb{s)/cr^f^{s), where (t°^ — iira^/Ss 



is the lowest-order cross section for e"^ 



and 



CTb is the total cross section for e~^e~ bb{-^), including 
bb states produced in initial state radiation (ISR) below 
the open beauty threshold, i.e. the T(IS'), T(2S'), and 
T{3S') resonances. 

The particles produced in the collisions are detected by 
the BaBAR detector, described elsewhere 0. Charged- 
particle tracking is provided by a five-layer silicon ver- 
tex tracker (SVT) and a 40-layer drift chamber (DCH). 
In addition to providing precise position information for 
tracking, the SVT and DCH also measure the specific 
ionization (dE/dx), which is used for particle identifi- 
cation of low-momentum charged particles. At higher 
momenta (p > 0.7 GeV/c) pions and kaons are identi- 
fied by Cherenkov radiation detected in a ring-imaging 
device (DIRC). The position and energy of neutral clus- 
ters (photons) are measured with an electromagnetic 
calorimeter (EMC) consisting of 6580 thallium-doped 
Csl crystals. These systems are mounted inside a 1.5- 
T solenoidal super-conducting magnet. Muon identifi- 



cation is provided by the magnetic flux return system 
instrumented with Resistive Plate Chambers and Lim- 
ited Streamer Tubes. The full detector is simulated, for 
background and efficiency studies, with a Monte Carlo 
program (MC) based on GEANT4 0. 

To measure Rb, we count the number of events pass- 
ing a selection that enriches the sample in events con- 
taining B mesons (Nh) and those passing an indepen- 
dent di-muon selection (N^) at each energy point and 
at a reference energy below the open beauty production 
threshold. Indicating with a prime the quantities at the 
reference energy, we write: 



Nhis) 



(Rbis) cr°^ (s) - cjisRis)) esis) 



^c^x(s) ex{s) + crisR{s)eisR{s) 



(1) 



X 



K 



N' 



E 

V, X 



'^'x ^'x + '^'iSR. ^'iSR 



(2) 

(3) 
(4) 



where es is the efficiency for open b production to satisfy 
the hadronic selection, X represents the different back- 
ground components described later, ai represents the 
cross-sections for the process i, ti the corresponding ef- 
ficiency, and C is the integrated luminosity collected at 
a given value of -^/s. Measurements of and N'^ are 
needed in order to normalize the hadronic rates to the 
collected luminosities. As reference we choose the sam- 
ple collected at v^=10.54GeV, about 40 MeV below the 
T(45') mass, taken during 2006-2007. Special mention 
is made of the ISR sample, the production of T{nS) 
{n = 1,2,3) mesons via initial state radiation: albeit 
part of the signal, this process can occur at the reference 
energy and has an efficiency and an energy dependence 
of the cross-section different from the open beauty pro- 
duction. 

Solving the system of equations one obtains: 
/iV,(,s) N'^ \ e^js)^^ 

= [Ju^) w/^^^'^ ) + 

where we defined: 



J2x Rx{s)ex{s) + Risr{s) eisRjs) 

Sjf ^'x^'x + ^'iSR ^'iSR 



, (6) 
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and Ri = (Ti/a°^^ for each process and = a-fj,^/a°^, 



as- 



sumed independent of ^/s. It should be noted that these 
equations assume that the background scales with the in- 
tegrated luminosity, i.e. that the machine background is 
negligible, and that the di-muon selection leaves a negli- 
gible level of background. 

We select the 6-enriched sample by requiring at least 
three tracks in the event, a total visible energy in the 
event greater than 4.5 GeV, and a vertex reconstructed 
from the observed charged tracks within 5 mm of the 
beam crossing point in the plane transverse to the beam 
axis and 6 cm along the beam axis. These quantities are 
computed using exclusively tracks in the fiducial volume 
of the DCH (i.e. forming an angle with the beam axis 
0.41 < 6 < 2.54 rad). A further rejection of the main 
backgrounds, e'^e~ — s- qq, q = u,d,s,c events ("contin- 
uum" events) and e+e^ £^£^ , £ = e,fi,T events, is ob- 
tained by means of a cut on the ratio of the second and ze- 
roth Fox- Wolfram moments ;8:] , R2 , calculated using only 
the charged tracks. After optimization of the statistical 
sensitivity, we require R2 < 0.2. Events that pass this se- 
lection at the reference energy comprise 91% continuum. 



2% two photon (e+e 



g+e 7*7* 



e+e-Xh), and 7% 



ISR (e^ 



T{nS)jiSR) events. 



To select di-muon events, we require that two tracks 
have an invariant mass greater than 7.5 GeV/c^; their 
angle with the beam axis in the center-of-mass frame, 
6cms, must satisfy cosOcms < 0.7485, and the two muons 
must be collinear to within 10°. To exploit the fact that 
muons are minimum ionizing particles, we require that at 
least one of them leaves a signal in the EMC, and neither 
deposits more than 1 GeV. 

In the following we describe the method used to derive 
the inputs to Eq. O and the corresponding errors, sepa- 
rating correlated and uncorrelated errors. The covariance 
matrix for the measurements of Rb at different energies is 

VtJ = [<^stat{Si) + CrlnciSi)] % "t-CTcorr (Si )<Tcor-r (Sj ) , whcrC 
O'stat 

, (Jcorr 

(si), and cr„„c(si) are the statistical, cor- 
related, and uncorrelated systematic error respectively, 
and Sij is the Kronecker delta. 

The efficiency for the di-muon selection is extracted 
from a sample of fully simulated MC events generated 
with KK2f [9j at several values of ^/s. Due to the change 
in boost this efficiency is found to change by 1.5% over 
the whole range and the MC statistics error we assign 
to the corresponding correction is 0.2%. The correlated 
uncertainty on the absolute scale of the efficiency is esti- 
mated to be 1% and to come primarily from uncertainties 
in the simulation of the trigger, of the quantities used 
in the selection and of the tracking efficiency. We also 
account for differences in the trigger configurations be- 
tween the scan data and the reference data taken during 
the year 2007 and estimate the efficiency on the reference 
data to be lower by (0.5 ± 0.2)%. The same generator is 
consistently used to extract = 1.48 ± 0.02, where this 
correlated error is due to the uncertainty on the cross- 



section. 

The efficiency for e'^e~ — > bb events is estimated by us- 
ing EvtCen [lO| as generator, separately for each possible 
two-body final state including B, Bg, and B* mesons, and 
at different values of y^. Because we ignore the relative 
composition in terms of final states at each energy we 
consider the largest and the smallest efficiencies among 
the allowed final states and take their mean value as the 
central value and half their difference as uncorrelated er- 
ror. The correlated error on the absolute scale of €b is 
estimated by varying the selection criteria and it is found 
to amount to 1.3%. 

The calculation of the double ratio k^^ requires the de- 
pendence on y/s of e^, which has already been discussed, 
and the cross-sections and efficiencies for the ISR and 
the background processes. 

The ISR cross-section is computed to second-order ac- 
cording to Ref . [ll| . The corresponding efficiency (e/s/j) 
is estimated with MC simulation to be 41% on average. 
The relative efficiency change across the scan, estimated 
to be ~ 5%, is used as a correlated uncertainty and it 
propagates to an error on Rt of at most 0.7%. 

The cross-section for two-photon events scales as the 
square of the logarithm of s, and the corresponding effi- 
ciency is considered to be flat. The product of the cross- 
section and the efficiency (a-y-ye-y^) before the R2 is fitted 
from the distribution of the direction of the missing mo- 
mentum and then multiplied by the R2 cut efficiency. We 
attribute 50% uncertainty to this estimate, leading to a 
relative correlated error of at most 0.2%. Finally, the 
product of the continuum cross-section and efficiency is 
computed by subtracting the ISR and two photon com- 
ponents from A''^ (see Eq. [2]). The continuum contribu- 
tion to R (Rcont) is assumed to be constant with -^i, 
while the corresponding efficiency (econt) was estimated 
on a sample of MC events generated with JETSET [l3| ■ 
No correction to account for the fact that the reference 
data were taken in a different data-taking period was 
found necessary. The relative change of econt over the 
whole scan range is estimated to be 3% and a 0.2% sys- 
tematic error due to MC statistics is assigned to it. We 
also find that the distribution of R2 is not perfectly re- 
produced by the MC. We therefore estimate the scahng 
of Ccont separately with and without the R2 < 0.2 re- 
quirement and take the difference among the results as a 
correlated systematic error. Its contribution depends on 
the value of Rb and it is at most 2%. 

To measure ^/s of each point we fit the distribution of 
the invariant mass of the two muons in the selected di- 
muon sample with a function made of a Gaussian with 
an exponential tail on the side below the peak mass. We 
then use the mean of the Gaussian as estimator of ^/s 
and we determine a bias of (20.9± 1.5) MeV for this quan- 
tity by comparing the T(3S') mass measured on the data 
taken during the ^ 100 pb~^ scan performed by PEP-II 
at the beginning of the last data-taking period with the 
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FIG. 1: (Left) Measured iij as a function of ^/s with the position of the opening thresholds of the e^e B^J^B 'J^ processes 
indicated by dotted lines. (Right) A zoom of the same plot with the result of the fit described in the text superimposed. The 
errors on data represent the statistical and the uncorrelated systematic errors added in quadrature. 



TABLE I; Contributions to the relative correlated systematic 
error on Rb- 



Contribution 


Relative error (%) 


HIJ, MC statistics 


0.2 


/i/i radiative corrections 


1.4 




1.3 




1.3 


^cont 


< 2.0 


^ISR 


< 0.7 




< 0.2 



resonant depolarization result Il3!|. We correct for this 
bias, that comes from the (strongly) non linear impact 
of the momentum resolution in the invariant mass, and 
verify on simulated events that it does not depend on y^. 



TABLE II: Fit results for the T(10860) and T(11020) res- 
onances resulting from the fit described in the text. The (j) 
phases are relative to the interfering continuum. The corre- 
sponding world averages 17] are also reported. 





T(10860) 


T(11020) 


mass (GeV) 
width (MeV) 
(f> (rad) 


10.876 ± 0.002 
43 ±4 
2.11 ±0.12 


10.996 ± 0.002 
37 ±3 
0.12 ±0.07 


PDG mass (GeV) 
PDG width (MeV) 


10.865 ± 0.008 
110 ± 13 


11.019 ±0.008 
79± 16 



The resulting measurements of as a function of ^/s 
are shown in Fig. [l] where the error bars represent the 
sum of the statistical and uncorrelated systematic errors 
and dotted lines show the different B meson production 
thresholds. The relative correlated systematic errors on 
Rh are summarized in Table [D The numerical results 



for each energy point, together with the estimated ISR 
cross section, can be found in Ref. It is important 
to stress that radiative corrections have not been applied 
since they would require an a-priori knowledge of the 
resonant region. The measured i?b therefore includes all 
final- or initial- state radiation processes. 

The large statistics and the small energy steps of this 
scan make it possible to observe clear structures corre- 
sponding to the opening of new thresholds: dips cor- 
responding to the B'-*^B* and BgB* openings and a 
plateau close to the B*B* one. It is also evident that the 
T(10860) and T(11020) behave differently above and be- 
low the corresponding peaks. Finally, the plateau above 
the T(11020) is clearly visible. 

We fit the following simple model to our data be- 
tween 10.80 and 11.20 GeV: a flat component represent- 
ing 66-continuum states not interfering with resonance 
decays, added incoherently to a second flat component 
interfering with two relativistic Breit Wigner resonances: 
a = \Anr\^ + \Ar + ^108606"^^°^™ BI^ (M10860 , Tioseo) + 
Aiio2oe''^"''^''W(Miio2o,riio2o)|', with BW{M,T) = 



l/[{s- AP)+iMT]. The results summarized in Table|TT] 
and Fig. [1] differ substantially from the PDG values. In 
particular the B*Bg and B*B* thresholds have a very 
large impact on the determination of the T( 10860) width. 

The number of states is, a priori, unknown as are their 
energy dependencies. Therefore, a proper coupled chan- 
nel approach 1^ including the effects of the various 
thresholds outlined earlier, would be likely to modify the 
results obtained from our simple fit. As an illustration 
of the systematic uncertainties arising from the assump- 
tions in our fit, a simple modification is to replace the flat 
nonresonant term by a threshold function at ^/s — 2m b- 
This leads to a larger width (74 ± 4 MeV) and a lower 
mass (10869 ± 2 MeV) for the T(10860). 
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In summary, we have performed an accurate measiirc- 
ment of Rh in fine grained center-of-mass energy steps 
and have shown that these measurements have the poten- 
tial to yield information on the bottomonium spectrum 
and possible exotic extensions. 
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